C§ axial form factor from bubble chamber experiments 



K. M. GraczyktQ D. Kielczewska^, P. Przewlocki*, and J. T. Sobczyk* 
^Institute of Theoretical Physics, Wroclaw University, pi. M. Borna 9, 50-204, Wroclaw, Poland 
* Institute of Experimental Physics, University of Warsaw, Hoza 69, 00-681 Warsaw, Poland 
* A. Soltan Institute for Nuclear Studies, Hoza 69, 00-681 Warsaw, Poland 
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agreement. For the dipole parametrization of C^(Q 2 ), we obtain C5 (0) = 1.19 ± 0.08, Ma = 
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current l7r° production cross section, important for the T2K neutrino oscillation experiment. 
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1. INTRODUCTION 

The problem of the correct evaluation of neutrino sin- 
gle pion production (SPP) cross sections is interesting 
by itself and also important for future neutrino oscil- 
lation experiments. Recent measurements suggest that 
SPP cross section may be 20%-25% larger than what was 
assumed in the past [1| . Understanding of nuclear effects 
in neutrino interactions in the 1 GeV energy region is 
still not satisfactory. In order to cope with the deficit 
of events in the low Q 2 region experimental groups have 
introduced effective quantities like k in charge current 
quasielastic scattering (CCQE) [2| or very large value of 
Fermi momentum in Pauli blocking 0] • These difficulties 
have the positive effect that neutrino cross section exper- 
imental results are more often presented in the form of 
raw data, which are only efficiency corrected. Such mea- 
surements are free from dependence on the models im- 
plemented in the Monte Carlo generators of events but 
extraction of free nucleon cross-sections requires good 
knowledge of final state interactions (FSI) effects. For 
example, in the recent MiniBooNE neutral current (NC) 
l7r° cross-section measurement on CH2 there were large 
carbon pion absorption and charge exchange nuclear ef- 
fects 

For this reason it makes sense to come back to old but 
relatively good statistics SPP data on deuteron obtained 
in bubble chamber experiments in the Argonne National 
Laboratory (ANL) la and in the Brookhaven National 
Laboratory (BNL) [6[ . The big advantage of these mea- 
surements is that deuteron nuclear effects are easier to 
control. It was claimed that SPP data from two exper- 
iments are in disagreement 0, 0], because of different 
total cross-sections and shapes of dajdQ 2 . Apparently 
it is true that the total cross-sections reported by the 
BNL experiment seem to be systematically larger than 



those obtained by the ANL. But it should be noted that 
in both experiments there are large normalization uncer- 
tainties of the measured cross-section due to imprecise 
knowledge of the neutrino flux. It is therefore desirable 
to perform a careful simultaneous reanalysis of two sets 
of data and this is the aim of our paper. We notice that 
during the last few years several theoretical models have 
been proposed 0, [a, 0, EH to describe the weak pion 
production off nucleon/nucleus. They were fine-tuned to 
the neutrino scattering data (usually only to the ANL 
data) not corrected to the effects which we consider in 
this paper. 

The details of our approach will be given in the next 
sections of the paper, and here we would like to outline 
the main points: 

1. We focus on the v^p — ► n~pir + reaction because for 
this process it seems to be reasonable to neglect a small 
nonresonant contribution. However, we notice that 
some of the current theoretical approaches include non- 
resonant background also for this channel (see [loj]). 

2. We will include deuteron nuclear effects applying the 
approach from Refs. 
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3. In our statistical analysis we use the \ 2 function with 
a contribution also coming from the overall normaliza- 
tion (flux) uncertainty, and this turns out to be impor- 
tant. 

Our main result is that ANL and BNL data are in 
quite a good agreement. We demonstrate this first by 
providing the value of \ 2 f° r our best fits but then also 
using more sophisticate methods of checking the self- 
consistency of the two independent sets of data. 

We present our results in the form of a fit for 
the C^(Q 2 ) form factor. We investigated several 
parametrizations with two options: either keeping C$(Q) 
as a free parameter or fixing its value using the argu- 
ment motivated by the partially conserved axial cur- 
rent (PCAC) hypothesis. In the most standard dipole 
parametrization with the Cd(0) kept as a free parameter 
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and with correction from deuteron effects included our 
best fit is: C A (0) = 1.19 ±0.08, M A = 0.94 ±0.03 GeV. 

We also present two applications for our results. The 
first one is motivated by a need to evaluate the NC 1tt° 
cross-section in the T2K experiment 14, [lij]. We will 
estimate the uncertainties of the NC 1tt° cross-section. 
In order to make the evaluation realistic, we will use two 
different Monte Carlo event generators with all the nu- 
clear effects included. The second application is a com- 
parison with recent measurements of the CCl7r + /CCQE 
ratio done by MiniBooNE on the CH 2 target . This 
measurement is particulary important because it is free 
from the flux normalization uncertainty. We obtain good 
agreement with the data. 

Our paper is organized as follows. Sec. 2 contains ba- 
sic definitions and notation, and in Sec. 3 we introduce 
various parametrizations of the axial form factors. In 
Sect. 4 we describe our approach to treat deuteron nu- 
clear effects. Sec. 5 is the most important, as wc present 
the methodology of our analysis: we define x 2 and for 
completeness we provide all the necessary information 
about the form in which ANL and BNL data are given 
and used in our numerical analysis. In Sec. 6 we show 
our results for various parameterizatons of the axial form 
factor. The complete set of results is given in Tables. U 
and [ill Sec. 6 contains also detailed discussion of our 
results as well as a formal demonstration that ANL and 
BNL data sets are self-consistent. In Sec. 7 applications 
of our results are given and last Sec. 8 contains our final 
remarks. 



2. ADLER-RARITA-SCHWINGER FORMALISM 

We analyze the experimental data for the charged cur- 
rent (CC) neutrino-proton reaction: 



where fc, k' , p and p' denote neutrino, muon, proton and 
A(1232) resonance four-momenta. The four-momentum 
transfer and its square are 



q = p' — p=k — k', Q 2 = -q 2 



and the hadronic invariant mass is 



W 2 = p' 2 = {p + qf. 



(2) 



(3) 



One way to describe the reaction ([I]) is to apply the 
Adler-Rarita-Schwinger formalism. The final hadronic 
state is a 3/2-spin resonance described as a Rarita- 
Schwinger field. The transition from the nucleon to 
A ++ is given as a matrix element of the weak hadronic 
current, which has the standard vector-axial structure 



J? C = JV+J£. (4) 
Both vector and axial parts are expressed in terms of 
several form factors. 



v{k) +p(p) -> n-(k') + A ++ (p') 



(1) 



Under general assumptions, the vector part can be 
expressed by means of three form factors, C^iQ 2 ), 
CY(Q 2 ), and CUQ 2 )-- 



(A++(p>)\jV\N(p))=V3* x ( P >) 



' r<v n v r<v \ /nv nv /-<v 

■ \ I + —i-r/ + il5_ B )a" -a x ^-t + -W + 



r 



75W(P)) 
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where M is the nucleon mass, ^(p') is the Rarita- The axial part depends on four form factors C^(Q 2 ), 
Schwinger field, and u{p) is the Dirac spinor. C^(Q 2 ), cf(Q 2 ) and C^(Q 2 ) [H|: 



(A++(p')\jt\N(p)) = yffixtf) 



A I w 3 i ~-4 I \ V A I "3 i "4 / I | A /-~f 



M 5 M 2 



■Ct 



u(p). 



(6) 
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The differential cross-section for the reaction ([1} reads current reaction G = G), and 9 C is the Cabibbo angle. 

The neutrino energy is denoted by E. 

Wfu, and L^v are hadronic and leptonic tensors defined 
where G = G cos 9 C , G is the Fermi constant (for a neutral as: 



»V - ^ E <A^y | Jg° IP) <A^y | J~ IP)' _ ^ , (8) 



L^iy — 8 ^fc ^fcjy + fc^fc j/ Qiu/k a k +~ i^^i/apk k J . (9) 



Ta(W / ) is the A width, for which we assume the P- 
wave (Z = 1) expression 

iw.fi^r 1 ^ a») 



q cm (M A )J W 



with 



Qcm (W) 



fW 2 +M 2 -ml 
\ 2W 



- M 2 (11) 



M A = 1232 MeV, ra„ = 139.57 MeV is the charged pion 
mass. 



Additionally, the amplitude A 1 a is computed like (fT2"j) . 
but as an input one needs to apply hadronic current ac- 
cording to the recipe from the standard model: the vector 
part is multiplied by (1 — sin 2 9w), @w is the Weinberg 
angle, and the axial part is left unchanged 

A' A = A A (CY -» (1 - sin 2 e w )CY, Cf -> Cf) . (19) 



3. FORM FACTORS 



3.1. Vector form factors 



2.1. Neutral current scattering 

We introduce the amplitude for the p — > A ++ weak 
transition: 



A A = A ( P -> A- 



• pTT 



(12) 



Then the amplitudes for the other two channels for CC 
neutrino-neutron scattering are expressed in terms of the 
appropriate Clebsch-Gordan coefficients: 



-^A A =A(? 



1 



A A =A (n -> A+ 



P7T 



(13) 



(14) 



A way to describe SPP in NC neutrino-nucleon scatter- 
ing was proposed in Ref. [lH . First of all, the charge cur- 
rent and neutral current amplitudes are related through 
the Clebsch-Gordan relations: 



A' (up - 


■* v/a+ - 


-* l/TT°p) — 


V2 

T ^ A 


(15) 


A'(vn - 


-+uA°- 


-> vif'n) = 


—-4 a 


(16) 


A\vp - 


+ i/A°- 


■> ^7r + ra) = 


3^ 


(17) 


A\vn - 


► v/a + - 


-> V7T~p) — 


lA>A. 


(18) 



In an older analysis (see e.g. [23j) an additional con- 
straint on the vector form factors was imposed: 

CY(Q 2 ) = 0, CY(Q 2 )=^CY(Q 2 ) (20) 

as motivated by the quark model relations [SU (6) sym- 
metry relation, for details see e.g. Ref. [13|]- The above 
relations describe the dominance of the magnetic ampli- 
tude Ml, while the electric amplitude E2 vanishes. With 
constraints given in (|20|) . the vector current is expressed 
by only one unknown function, which can be extracted 
from the electroproduction data. 

In the original analysis of the ANL and BNL data two 
parametrizations of the (Q 2 ) were considered. The 
first one proposed by Dufner and Tsai j2f| (applied to 
ANL data analysis [26j) reads 



C 3 y (Q 2 ) = 2.05v/TT9Qe- 315 «, Q=y/Q 2 , 

and the other one is of a simple dipole form, 

2.05 



c v AQ z ) = 



1 



Q 2 



(21) 



(22) 



0.54 GeV 2 , 

In our analysis we will use instead the recent experi- 
mental fits proposed in [271 ]: 

CY(Q 2 ) = 2.13(1 + St) G D (Q 2 ), (23) 



4M 2 
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CY(Q 2 ) = —1.5 i [ I 
CY{Q 2 ) = 0.48 ( I 



Q 2 



AM 2 

Q 2 



0.776M 2 



Gd(Q 2 ), (24) 



Gd(Q 2 ), (25) 



where 



and Afv = 0.84 GeV. 



(26) 

It is clear that the use of different vector form factor 
has an impact on the axial contribution as obtained from 
a fitting procedure to the neutrino scattering data. We 
will return to this point in the discussion. 



In the models presented in Refs. [13, HI, E3| , the value 
of C A (0) is smaller than the one obtained from formula 
(f2"9"| . Therefore it seems reasonable to use the experi- 
mental data in order to evaluate its value. In this paper 
we perform a comprehensive fit to the existing data on el- 
ementary interactions from ANL and BNL experiments. 
The value of CHO) was already studied experimentally 
by Barish et al. [26| with the conclusion that C A (0) = 1.2 
is in agreement with the data with an accuracy of 20%. 

In recent years the functional form of C ^ ( Q 2 ) has been 
also a subject of intensive research 0, Efl E3, H, [H, 0, 
H, EH, [H, A possible parametrization of C A (Q 2 ) is 
based on a comparison of the Rarita-Schwinger formalism 
with the predictions of the Adler model [23j : 



3.2. 



Axial form factors 



C A (Q 2 ) = C A (0) 1 + 



aiQ 2 

h + Q 2 



l + (30) 



The main contribution to the axial current comes from 
C A (Q 2 ), which is an analog of (Q 2 ) in the vector 
current - it describes the (Ml)" 4 amplitude. Similarly, 
C A (Q 2 ) is an analog of C§ (Q 2 ) because it contributes to 
the (E2) A amplitude In the SU(6) symmetry limit 
(E2) A = 0, C A {Q 2 ) = C A (Q 2 ) = and only C 4 A (Q 2 ) 
and C A (Q 2 ) are nonvanishing. 

The SU(6) symmetry relations are only approximate, 
and they become broken by color hyperfine interaction. 
A lot of effort was done to compute all axial form factors 
directly from the quark models [28[ but the results do 
not reproduce the data sufficiently well. 

In the phenomenological analysis there is no enough 
experimental data to extract all the axial form factors 
separately. Therefore some extra constraints must be 
imposed: 

(i) typically one sets C A (Q 2 ) = 0; 

(ii) the Adler model [l8| suggests 



with values: 



C 3 » 



«4 



a-3 = 
1.2, 

a 5 = 



6 3 = 0, Cf(0) = -C 5 A (0)/4, 



-1.21, b 4 = b 5 = 2. (31) 



ofC A 



In our analysis the first ansatz for the functional form 
(Q 2 ) will be as above. The fitted parameters will 
be either only Ma or both Ma and C A (0). 

In the second fit we consider the simplest possible func- 
tional form, namely the dipole parametrization: 



C A (Q 2 ) = 



cm 



M\ 



C A (0) = 1.15. 



(32) 



Again, we will fit either only Ma or both Ma and C A (0). 

Let us emphasize that in both cases when C A (0) was 
treated as a free parameter we kept the constraint ([2T)l . 



Ci(Q 2 



-C£(Q 2 )/4; 



(27) 



NEUTRINO-DEUTERON SCATTERING 



(iii) the PCAC hypothesis implies 

M 2 
m%. + Q z 

with m-n- being the pion mass. 



C A (Q 2 ) 



(28) 



C5 (0) can be evaluated from the off-dia gon al 
Goldberger-Treiman relation (for a review see e.g. [29(). 
The updated value is reported in Ref. [28j : 



V6M 



= 1.15 ±0.01, 



(29) 



where the values g^NAiq 2 — m 2 .) — 28.6 ± 0.3 and f v 
92.4 MeV were used. 



In both the ANL and BNL experiments for most of the 
exposition the detectors were filled with deuteron. The 
nuclear effects for the A ++ production were discussed in 
Ref. The final state interactions can be neglected 

and the effect comes from the fact that the target proton 
is bound. 

In [l2T ] the deuteron wave functions were considered to 
be obtained in three different nuclear models: Hulthcn 



34( , Paris [35j , and Bonn [36J . The results of [12[ indicate 
that the nuclear effects for A ++ production are larger 
than for the quasielastic scattering, in which case the 
main modification is a reduction of the cross-section in 
the region of small Q 2 < 0.05 GeV 2 due to Pauli blocking 
while nuclear effects are negligible for Q 2 > 0.15 GeV 2 

37 1 . According to the results of [l2j in the case of A ++ 
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Q 2 [GeV 2 ] Q 2 [GeV 2 ] 

FIG. 1: Differential cross-sections for v + d — > |ti _ + n + p + ir + scattering computed for the ANL beam. The solid/dashed 
lines denotes the best fit obtained with C^(0) fitted/fixed (with value 1.15). In the left panel the cross-sections computed with 
the dipole functional form of C£(Q 2 ) ^ are shown, while in the right panel results obtained with Adler parametrizations 
(|30|) are presented. The black squares denote the experimental data [jj, while the black triangles denote the experimental data 
multiplied by factor of 1.08 (left panel) and 1.11 (right panel) - the renormalization factors were obtained from the global fits 
(see Tables U and [TTJ . The theoretical cross-sections are modified to include deuteron nuclear effects. The cut W < 1.4 GeV on 
the hadronic invariant mass was imposed. 



excitation the nuclear effects are slowly varying with Q 2 
and they reduce the differential cross-section in Q 2 by ~ 
10% in the case of the ANL beam and ~ 5% for neutrinos 
of energy 1.6 GeV (a typical value for the BNL beam). 

In our analysis we therefore assume that the neutrino- 
deuteron A ++ excitation differential cross-section in Q 2 
gets modified by a function R(Q 2 ) with respect to the 
neutrino-proton cross-section: 



da 



deuteron 



(33) 



free target 



In general R{Q 2 ) is a function of the neutrino energy E, 
but for the ANL experiment we use the result of Ref. [l2j 
(Fig. 5), where the flux averaged differential cross-section 
is plotted with and without deuteron effects. Based on 
that we extracted: 



Ranl{Q 2 



{da{vd -> fi-nA++)/dQ 2 ) 



deuteron 



(da(vp^^A++)/dQ 2 ) freetarget 



(34) 



In the case of the BNL experiment we use R(Q 2 ) eval- 
uated at E = 1.6 GeV. 



There is an additional ambiguity in applying the re- 
sults from [T2| because the calculated reduction of the 
cross-section depends on the model of the deuterium po- 
tential. 

We considered all three nucleon-nucleon potentials 
mentioned above. However, in this paper only results 
obtained with the Paris potential are presented. The 
Bonn potential gives rise to very similar results. With 
the Hulthen potential, the results are still comparable, 
but the impact of nuclear effects is smaller. 

In [38| early ANL data on mostly hydrogen targets i.e. 
without nuclear effects, are presented. There are alto- 
gether 153 such events out of which it is estimated that 
105 are on the hydrogen target. This sample of events 
was used in the analysis of Schreiner and Von Hippel [23| . 



5. REANALYSIS OF THE BUBBLE CHAMBER 
DEUTERIUM DATA 

We look for a simultaneous fit to both ANL and BNL 
data by applying the x 2 method. The best fit is obtained 
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FIG. 2: Distribution of the number of events of the BNL experiment for v-\-d — > fi~ + n + p + -K + scattering. The solid/dashed 
lines denotes the best fit obtained with C5 (0) fitted/fixed (with a value of 1.15). In the left panel the cross-sections computed 
with the dipole functional form of C^(Q 2 ) (|32l) are shown, while in the right panel results obtained with Adler parametrizations 
(|30l) are presented. The black dots denote the experimental data [g|. The theoretical cross-sections are modified to include 
deuteron nuclear effects. The cut W < 1.4 GeV on the hadronic invariant mass was imposed. For higher Q 2 (above 1.15 GeV) 
the data bins were joined, to get better statistics. The new Q 2 bins have the following widths: 0.1 for Q 2 £ (1-15, 1.55), 
0.2 GeV 2 for Q 2 g (1.55, 1.95) GeV 2 , 0.4 GeV 2 for Q 2 e (1.95, 2.35) GeV 2 , and 0.65 GeV 2 for Q 2 <= (2.35, 3) GeV 2 . 



by minimizing a function: 



X' = Xanl + Xbnl> (35) 

where Xanl an< ^ X 2 bnl are defined for each data set sep- 
arately. 



In both cases the % 2 is given by the standard formula 
with an additional quadratic term which comes from the 
total systematic uncertainty for the flux [4o| : 



* 2 = E 



AN, 



I a tot-th 



N, 



J^th 



(36) 



or equivalently by 



" f„diff fn2 



p-l 



(37) 



with: 



'tot — th 



17 N th ' 

tot — exp 



(38) 



N { and N th i are experimental results and theoretical 
predictions for the number of events in the i-th Q 2 -bin, 



AA^ is the experimental result uncertainty (it is the sum 
of statistical and systematic uncorrelated contributions) , 
a tot-exp an d o'tot-th are the experimental and theoreti- 
cal flux averaged cross-sections measured and calculated 
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Adler 
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15 (from PCAC) 
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FIG. 3: (7^(0) fits and la error bounds obtained by assuming 
dipole [Eq. (|32[) ] and Adler [Eq. (|30| parametrizations of axial 
form factors. The solid perpendicular line denotes the PCAC 
value quoted in f2f| together with la error bounds. 



with the same cuts, and finally 



N - 



N 



th 



E^ 



th,j- 



(39) 



In both above cases the form factor parameters are 
fitted but in the first case [Eq. (|36)) ] the overall number 
of events N t h is fitted, while in the second case (Eq. ([37]) ) 
the fitting is applied to p. In the presentation of our 
results, in both cases the final results are presented by 
giving the values of p. 



5.1. ANL data 

The ANL data are given in the form of normalized 
da/dQ 2 @] with a cut W < 1.4 GeV. Information about 



the ANL neutrino flux is provided in Ref. [39j (see Fig. 
8). The uncertainties in the differential cross-section 
(both statistical and uncorrelated systematic) are given 
in the original paper. 

In the case of ANL experiment the data consists of 1115 
(871) corrected (raw) events. We use the flux averaged 



da 



dQ^ Ql) 



2^ _ a ANL,diff 



(.//' 



(40) 



data points (all are taken from Radecky et al. (H) with 
uncertainties A (^j^ NL ' dl ff^j . All the numbers are given 

in table IV of Ref. [5j. Cuts are imposed: for neutrino 
energy E 6 (0.5,6) GeV and for the hadronic invariant 
mass W < 1.4 GeV. The data covers the range in Q 2 
from Q 2 = 0.01 GeV 2 to Q 2 = 1 GeV 2 

The experimentally measured (flux averaged and with 
the cut TV < 1.4 GeV) cross-section is calculated to be: 



ANL 



E 



ANL,diff^Q2 



i=l 



0.31 x 10~ 38 cm 2 , 
(41) 



where AQ~ are the bin widths, Uanl = 9. 

The theoretical formula for the differential cross- 
section in a given Q 2 bin is the following: 



ANL,diff 
7 th 



(Ql) 



^A 



NL 



AQ 2 



dQ 2 

??-AQf/2 ./£',,. 



dE 



1.4 GeV 



dW<£ ANL (E)a th (E, Q 2 ,W), 



M+m„ 



r 



(42) 



with 



NL 



dE1> 



ANL 



(43) 



Emm = 0.5 GeV and E max = 6.0 GeV. 

The systematic uncertainty of the total cross-section 
resulting from an imprecise determination of the neutrino 
flux is quoted in Ref. [5] to be 15% for E € (0.5, 1.5) GeV 
and 25% for E > 1.5 GeV. The flux was calculated on the 
basis of pion production cross-sections measured by Cho 
et al. [56j for the same proton beam in a separate experi- 
ment. The calculation is described in Ref. [39( where the 
neutrino flux computed based on pion production was 



compared with the flux derived from measurements of 
QE interactions. The conclusion was that the latter flux 
was smaller by 21%. Therefore in our discussion, we as- 
sume the average overall normalization uncertainty to be 
20%. 



r ANL = 0.20. 



(44) 



In the case of ANL data it is natural to use formula ([3T)) 
for x 2 ■ 
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Ma (Co\i r \ 

1V1A V ) 


°5 yJ) 


PAN L 


Pbnl 


v 2 / WDP 
X /IV UP 


VjrOr 


dipole, only Ma, free target 


0.95 ± 0.04 




1.15 ±0.06 


0.98 ±0.03 


25.5/28 


0.60 


dipole, only Ma, deuteron 


U.94 ± U.04 




1.U4 ± U.OO 


U.9f ± U.03 


Z4.0/20 


U.6o 


dipole, Ma and C^(0), free target 


0.95 ± 0.04 


1.14 ±0.08 


1.15 ±0.11 


0.98 ±0.03 


25.5/27 


0.54 


dipole, Ma and C$(0), deuteron 


0.94 ±0.03 


1.19 ±0.08 


1.08 ±0.10 


0.98 ±0.03 


24.3/27 


0.60 



TABLE I: The results obtained for fitting the dipole parametrization [Eq. (|32f) ] of the axial form factor. 





M A (GeV) 


C 5 A (0) 


PANL 


Pbnl 


X 2 /NDF 


GoF 


Adler, only Ma, free target 


1.31 ±0.06 




1.23 ±0.06 


0.98 ±0.03 


27.1/28 


0.50 


Adler, only Ma, deuteron 


1.29 ±0.07 




1.11 ±0.06 


0.98 ±0.03 


24.8/28 


0.64 


Adler, M A and C^(0), free target 


1.31 ±0.07 


1.1 ±0.08 


1.18 ±0.1 


0.98 ±0.03 


26.7/27 


0.48 


Adler, M A and C^(0), deuteron 


1.29 ±0.07 


1.14 ±0.08 


1.11 ±0.10 


0.98 ±0.03 


24.8/27 


0.58 



TABLE II: The results obtained for fitting the Adler parametrization [Eq. (|30[l ] of the axial form factor. 



5.2. BNL data 

The BNL data are given in the form of the distribution 
of events in Q 2 [4l| with a cut W < 1.4 GeV. Information 
about the flux is given in [42| (Fig. 7). Reconstruction of 
the neutrino spectrum is also presented in a later reanal- 
ysis done by Furuno et al. [44]. In our analysis we use 
the flux from this latest paper (Fig. 1, right panel). The 
overall normalization (cross-section) is also provided but 
without cuts on the hadronic mass [4l|. The beams in 
the ANL and BNL experiments are quite different so that 
the results are expected to give independent information 
on da/dQ 2 . 

The BNL data sample consists of 1803 (1610) corrected 
(raw) events [4l[ . The statistics is better than in the ANL 
experiment by a factor of ~ 50% . In the case of the BNL 
experiment we consider dN/dQ 2 taken from Fig. 10 of 
Ref. [il| with only statistical errors ANi = \fWl. The 
events are collected under the condition W < 1.4 GeV 
and E e (0.5,6) GeV for Q 2 e (0,3) GeV 2 . For Q 2 > 
1.2 GeV 2 we decided to combine some bins in order to 




get better statistics (as illustrated in Fig. We can do 
this because in the analysis only the statistical errors of 
the experimental points are taken into account. 

The neutrino flux presented in the BNL papers 
was determined from measurements of QE events [42 
The QE axial mass derived only from the shape of Q 2 dis- 
tribution with precision of 6% was used to calculate the 
total cross-section, and the observed rate of QE events 
allowed one to then calculate the neutrino flux with an 
uncertainty of less than 10%. The axial mass was later 
recalculated in Ref. [57} with new electro-magnetic form- 
factors and the updated value was found to be only 2% 
smaller. The calculation of the neutrino flux coming 
from the same proton beam of 29 GeV/c from Alter- 
nating Gradient Synchrotron (AGS) was presented in a 
paper by Ahrens et al. [43| and compared with the flux 
derived from QE events in a larger scintillator detector, 
which could measure the beam profile, contrary to the 
deuterium bubble chamber. It has been found that both 
fluxes differ by 10% and the spectrum shapes agree very 
well. 
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Thus in our discussion below, we assume the normal- 
ization uncertainty of 10% and in the main analysis we 
take 

r BNL = 0.10 (45) 

(In Sec. 6.2 we will consider the impact of different 
normalization uncertainties on the final results.) To fit 
the normalization parameter for the BNL data we use 
the total cross-section data published in 41]. For the 
total cross-section data the hadronic invariant mass is 
unconstrained, and the neutrino energy range is £ £ 
(0.5,3) GeV. 

In its analysis, the BNL collaboration used only the 
events satisfying Q 2 > 0.1 GeV 2 . Two different justifica- 
tions for that can be found. In Q it is written that the 
efficiency of reconstructing events for Q 2 < 0.1 GeV 2 is 
very low. This point is later investigated in detail and 
the dependence of the reconstructed parameters (axial 
mass) on Q 2 min is shown. The value of the fit to the axial 
mass is rather stable for Q 2 nin > 0.06 GeV 2 . In [6j one 

I 



6. RESULTS AND DISCUSSION 

We analyze simultaneously both the ANL and BNL 
data. The number of degrees of freedom is 

NDF = n ANL + riBNL - n par - 2, (51) 

where n par is the number of parameters in the analytical 
expression for the C^(Q 2 ). The extra factor of 2 comes 
from two renormalization constants: panl and pbnl- 

As explained in the Introduction, two parametrization 
of C^(Q 2 ) are studied: 

(i) dipole [Eq. (32)] , 

(ii) Adler [Eq. l[50])]. 

We start the numerical analysis with discussion of the 
simplest parametrization of the axial form factor, namely 



can also find a justification that nuclear effects play an 
important role in the Q 2 < 0.05 GeV 2 region. 

In the reanalysis of the BNL data done in Ref. [44| it 
is stressed that -^jj data are presented with no nuclear 
corrections. 

A natural question arises as to why the ANL collabora- 
tion did not introduce a similar Q 2 cut on their data. A 
possible explanation is that in their case a Q 2 > 0.1 GeV 2 
cut would eliminate as many as ~ 15% of the events. This 
is because the ANL neutrino beam is of lower energy. 

Since the BNL data are given in a different form than 
in the case of ANL it is natural to use the expression 

The theoretical numbers of events in a given Q 2 bin 
are computed with the following formula: 

N BNL _ N th ^i a th W.i l A a.\ 

^th.i — ^ BNL BNL ' ^ °> 
a tot-th 

with 



(47) 
(48) 



(49) 
(50) 



the dipole one f|32[) . The results are summarized in Ta- 
ble U There are four different fits: with C^(0) = 1.15 or 
with C5 (0) treated as a free parameter and in both cases 
with and without deuteron corrections. All fits have ac- 
ceptable goodness. 

In the case what we consider the most reliable, i.e. 
with C§ (0) treated as a free parameter and with deuteron 
effects included we obtain 

Ma = 0.94 ± 0.03 GeV, (7^(0) = 1.19 ± 0.08. (52) 

It is interesting to see that in all four cases we obtained 
similar values for the axial mass. 

The inclusion of the deuteron effects does not affect 
much the theoretical parameters: Ma and C^(0). This 
is a consequence of the fact that in wide Q 2 range the 
deuteron effects mainly change the normalization of the 



-, i r Q 2 +AQ 2 /2 pE max r 1.4 GeV 

af h NL4in {Ql) = -r^ dQ 2 $ BNL (E)dE dWa th (E,Q 2 ,W), 

*BNL &Qi jQi-AQ 2 /2 JE min JM+m„ 

i r fE max r 1.4 GeV 

BNL _ 1 / in.2 I jp^, /en / jti/_ _Cc r\2 



'tot-th 



1 r rt^max LreV 

dQ 2 dE<5>BN L {E) dWa th (E,Q 2 ,W). 

*BNL J JE mi „ JM+m^ 



' E min J M+m„ 

In the case of BNL data the total cross-section is given without any cut on W and we need: 

"E„ 



°bnl = ~fr~~ — / " <5>BN L {E)<J exp (E)dE = 0Mx HT 38 cm 2 , 

^ BNL JE min 

"BNL = I dQ 2 ( dE<S>ANL(E)<J t h(E,Q 2 ). 

^ BNL J Jb„„ 



r 
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FIG. 4: Uncertainties of the fitted parameters. On the left the la error contour in the Ma and C^(0) plane is plotted. 
The black squares and triangles denote the points on the contour with maximal and minimal errors for total cross-section 
for v + p — » n~ + p + 7r + at different values of neutrino energy E. The obtained points are then mapped in the form of 
(C5 A (0)) ra i„,moi and (MA)min,max dependence on E (middle and right figures). 



da/dQ 2 reducing it by about 10%, and its impact be- 
comes compensated by the ANL renormalization param- 
eter panl- The normalization of the BNL data is not 
affected by the nuclear correction and the goodness of fit 
with the deuteron effects included becomes better. 

A similar analysis has been repeated for the Adler 
parametrization of C^(Q 2 ). The results are presented 
in Table ITT1 As before, there are four different fits: with 
C^(0) = 1.15 or C^(O) treated as a free parameter and 
in both cases with and without deuteron corrections. All 
fits have acceptable goodness. Again, the inclusion of 
deuteron structure correction affects mainly the normal- 
ization of the ANL data and leads to modification of 
C 5 » by ~ 3%. 

In the case that we consider the most reliable i.e. with 
C^(0) treated as a free parameter and with deuteron 
effects included we obtain: 

M A = 1.29 ±0.07 GeV, C^(0) = 1.14 ± 0.08. (53) 

In Fig. Q]the plots of da / 'dQ 2 {v^d fiA ++ ) are com- 
pared with the ANL data. The experimental points are 



re-normalized by factor panl according to the logic of 
our analysis. The theoretical predictions are multiplied 
by R{Q 2 ) in order to account for nuclear effects. 

In Fig. [2] the same is done for the BNL data. Above 
1.15 GeV 2 some bins are joined together, as explained in 
the introduction. 

We analyzed two different parametrizations of the 
C$(Q 2 ) form factor, and both led to the values of C$(0) 
which are compatible with PCAC result (see Fig. [3]). 

The quark models predict the direct relationship be- 
tween nucleon and resonance axial form factors [23| . For 
example in the Rein-Sehgal model [54l ] (an approach com- 
monly applied to neutrino data analysis) the nucleon and 
P33(1232) axial form factors have the same Q 2 depen- 
dence (see e.g. [32j), which, in practice, means the same 
axial mass. From that point of view it is interesting to 
notice that the values obtained by us for the axial mass 
of the dipole form factor are very similar to axial mass 
parameters extracted from the quasielastic bubble cham- 
ber data in the BNL {M% E = 1.07 ± 0.06 GeV) and 
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FIG. 5: Total cross-section for u-\-p — > p~ + p + 7r + . In the left panel the ANL data with the cut W = 1.4 are shown (black 
squares), while the right panel presents the BNL data |4ll ] (without cut in W) - black triangles. The overall normalization 
error is not plotted. The best fit curves were obtained with a corresponding cut in W . The theoretical curves were obtained 
with dipole parametrization Eq. (|32p with Ma = 0.94 GeV and C^(0) = 1.19. The shaded areas denote the la uncertainties 
of the best fit. The theoretical curves are not modified by the deuteron correction effect. 



in the ANL 5£ 



(M% E = 1.00 ±0.05 GeV) experiments. 



We examined the impact of a choice of vector form 
factors ([25]) on the axial mass Ma and on C^(0). Wc 
compared the dipole fits obtained by assuming two 
different parametrizations. The following numbers were 
obtained: 



Vector FF 


M A (GeV) 


cf(o) 


x 2 


PANL 


Pbnl 


Equation ([22]) 


0.93 


1.24 


30.6 


1.05 


0.97 


Equation {25} 


0.93 


1.19 


24.3 


1.08 


0.98 



We see that the axial mass does not depend much on 
the vector form factors which cannot be said about 
C^(0). For the vector form- factor (|22|) the higher value 
of C^(0) is compensated by the smaller panl but the 
shape of the obtained reconstruction of the differential 
cross-section is worse. 

Eventually, we fit the dipole form-factor to the BNL 
data, but we vary the Q 2 ut value. 



Qlut (GeV 2 ) 


M A (GeV) 


cm 


PBNL 


x 2 


GoF 


0.00 


1.01 


0.94 


0.96 


40.6 


0.05 


0.05 


0.96 


1.14 


0.98 


26.8 


0.47 


0.10 


0.93 


1.21 


0.98 


23.3 


0.67 



The dependence of obtained parameters on Q 2 ut is sig- 
As Q 2 cut becomes smaller the quality of the fit 



nificant 
becomes worse 



cussed here and the obtained values for Q 2 cut = 0.1 GeV 2 
can be recognized as the black triangle in Fig. 



6.1. Uncertainties of the fit and statistical 
consistency of ANL and BNL data 

The uncertainties of our fits and their impact on the 
uncertainties of cross-sections are analyzed by applying 
the covariance matrix algorithm [45| . The example of the 
covariance matrix for the fit with dipole parametrization 
(with deuteron correction) is given below: 



V = 10" 



( 



V; 



1.27 
-0.23 

1.19 
.2 • 10- 



-0.23 
7.41 
7.28 
0.44 



1.19 
7.28 
0.46 
0.54 



.2- 10- 2 > 
0.44 
0.54 
0.76 



(54) 



The correlation parameters are 
- -0.07, p 13 

P2i = 0.19, p 34 

of the 



P12 

pu = 0.08, 



= 0.33, P 23 = 0.83, 

0.19, (55) 



parameters is: 



We note that only BNL data are dis- 



where the sequence 

(M A , C£ (0) , PANL , PB NL ) ■ 

In Fig. [4] the la error contour in the M A and C^(O) 
plane is plotted. The maximal and minimal deviations 
of the cross-section values from the best fit, due to la- 
error are expected to lay somewhere on the error ellipse. 
For different values of the neutrino energy the points 
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0.85 0.9 0.95 1 1.05 

M A [GeV] 



(especially for the ANL), the second effect enlarges the 
data error, remembering that, one can say that the fit 
is compatible with the total cross-section data on the la 
level except from one ANL point which is in clear dis- 
agreement with all other measurements. 

In the case of the global analysis of the data obtained 
from independent measurements, the classical Pearson \ 2 
tests need not work efficiently. In Ref. [46j the statistical 
test, called parameter goodness of fit, devoted to verifica- 
tion of the statistical compatibility of different measure- 
ments is proposed. In the Appendix a short summary of 
this method is presented. The parameter goodness of fit 
test measures how far from the global minimum (of Xtot) 
the minima of separate data sets are. Applying this to 
our analysis we obtain [see Eq. IA2j 

X 2 = 0.2, NDF C = 2, (56) 

and the parameter goodness of fit PGoF = 0.9 [see Eq. 
([Aljl], This means that ANL and BNL data arc fully 
consistent. 

To illustrate our results we present in Fig. [6]the global 
minimum (denoted by an open circle) with la error el- 
lipse and the minima obtained separately for the ANL 
and BNL data, which are shown by black triangle and 
square (with la error bars) respectively. All points lie in 
the close neighborhood. 

6.2. Impact of flux uncertainties 



FIG. 6: The la (solid line), 2a (dotted line), and 3a (dashed 
line) error contours for the C^(0) and Ma fit parameters [see 
Eq. [32]. The best fit is denoted by the open circle. The 
best fit obtained for only ANL/BNL data is represented by 
black square/triangle. The open triangle denotes the best 
fit obtained with panl = Pbnl = 15% (the la contour is 
plotted with a dashed line). The open diamond denotes the 
best fit obtained with panl = Pbnl = 20% (the la contour 
is plotted with the dashed-dotted line). 



are placed in different region of the ellipse. We will use 
these results for the estimation of the uncertainties of the 
cross-sections obtained with theoretical predictions and 
the Monte Carlo simulations. 

In Fig. [5] the total cross-section for the v + p — » 
fi + +p+ir + is shown together with la uncertainties which 
are shown to be of the order of 10 %. The theoretical 
computation (with the dipole axial form factor) is com- 
pared with total cross-sections measured at the ANL [B} 
(left panel) and BNL [4l[ (right panel) experiments. We 
present ANL and BNL data on separate plots since the 
first are obtained with W cu t = 1.4 and the latter with no 
cut in W. The theoretical results are not corrected by 
the deuteron effect. We do not plot the overall normaliza- 
tion error. The first effect lowers the total cross-sections 



One can question whether the 10% normalization un- 
certainty for the BNL neutrino energy spectrum is not 
too optimistic. 

We checked how much our results depend on this as- 
sumption and we repeated our computations assuming 
both uncertainties to be first 15% and then 20%. For the 
dipole fit with nuclear effects included, we obtained the 
following: 

(a) 15% uncertainties: 



M A 


= 0.93 ±0.04 GeV, 


c£(o) 


= 1.16 ±0.10, 


PANL 


= 1.03 ±0.11, 


PBNL 


= 0.98 ±0.03, 


X 2 /NDF 


= 24.2/28, C.L. = 0.67 


(b) 20% uncertainties 




M A = 


0.93 ±0.04 GeV, 


C 5 » = 


1.15±0.12, 


PANL = 


1.02 ±0.14, 


PBNL = 


0.98 ±0.03, 



X 2 /NDF = 24.0/28 = 0.86, C.L. = 0.68. (58) 
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In Fig. [6] we added la curves for these two cases. 

We see that the considered modifications of the nor- 
malization errors only very weakly affect the obtained fit 
for the axial mass. The new central value of Crf(O) is in 
agreement with the one obtained before but the uncer- 
tainty becomes much larger. 

A general observation is that the increased normaliza- 
tion uncertainty enlarges the uncertainties of the form- 
factor parameters, and consequently, the cross-section 
uncertainties. However, we stress that rescaling of the 
normalization errors does not destroy the statistical con- 
sistency of the ANL and BNL data. 



7. APPLICATIONS 
7.1. Evaluation of NC 1ty° cross-section uncertainty 



It is interesting to investigate what the la error con- 
tours look like for different it production channels in 
terms of cross-sections. In particular, one might study 
two channels that are most important in modern neu- 
trino experiments, i.e. (i) + p — > /i~ + it + + p, 
+ n — > + 7r + + n (CC 7r + production) and (ii) 

v n + P ~~ * v n + n ° + Pi v n + n ~~ > v n + 71-0 + n ( N C n ° Pro- 
duction) . The 7T° production is the main source of back- 
ground in water Cherenkov far detectors of long baseline 
neutrino experiments searching for v e appearance, like 
T2K 0, [nl. The NC tt° events arc however difficult 
to study exclusively; one can try to study them by mea- 
suring 7r + production and extrapolating the results to 7r° 
production. 

In this section we will use fit results obtained previ- 
ously to estimate cross-section uncertainty for 7r produc- 
tion channels in the context of T2K experiment. For this 
purpose two software packages for simulation of neutrino 
interactions will be used. 

The main simulation package used in this analysis was 
the NuWro Monte Carlo event generator 47| . For reso- 
nant In production this generator uses the Adler-Rarita- 
Schwinger formalism for the A excitation. The non- 
resonant part is described by fraction of deep inelastic 
scattering (DIS) contribution, applying the algorithm de- 
scribed in [481 ] . NuWro is a generator in which many pa- 
rameters can be specified manually, including axial form- 
factor parameters, which are of interest in this analysis. 
Nuclear effects in oxygen (NEO) have recently been im- 
plemented in NuWro and were used in this study. 

As a reference, the Nuance neutrino generator [4!| was 
also used. Nuance is a widely used tool, tested in ex- 
periments with water Cherenkov detectors like K2K and 
Super-Kamiokande (SK). This appears to be consistent 
with the measurements of ir° production in the 1KT near 
detector of K2K [H(J [5l[ , as well as atmospheric neutri- 
nos in SK [52]. Its implementation of nuclear effects in 



oxygen can therefore be considered trustworthy (see Ref. 
[53j). Resonant pion production in this generator is cal- 
culated according to the Rein-Sehgal model [BH ]. 

We decided to examine the fit that used dipolc 
parametrization of axial form factors in which two pa- 
rameters Ma and C$(0) were fitted (see Sec. [3]). The 
error ellipse, presented in Fig. [4j is calculated for A ++ . 
However, we can scan this ellipse, calculating NC ir° pro- 
duction cross-section for each point on it, and find the 
minimum and maximum values (which correspond to the 
la range). By doing this for a broad range of incident 
neutrino energy we can obtain Ma and (0) parameters 
corresponding to the lower and upper la bounds at each 
energy and then use them in NuWro simulations, which 
will allow us to calculate cross-sections for the channels 
of interest. The NuWro simulation package uses exactly 
the same form factors as the ones that were used in the 
fit. 

All simulation samples in this work were created using 
water as a target - the most suitable material when simu- 
lating interactions in Super-Kamiokande, the far detector 
of the T2K experiment. Only muon neutrino interactions 
were taken into account as they dominate the T2K beam 

m. 

Figure [7] shows how cross-section uncertainty depends 
on energy of incoming neutrino. Two sets of points illus- 
trate la bounds on cross-sections (cf. Fig. obtained 
using NuWro; others are Nuance results. In order to show 
how the inclusion of NEO modifies the cross-sections we 
present separately results with nuclear effects turned off. 
All comparisons were done in the two 7r production chan- 
nels described earlier. Lower plots in Fig. [7] show all 
cross-sections divided by Nuance results without NEO. 
It is seen that the differences are most notable in the low 
energy region. In particular at 1 GeV the uncertainty for 
7T° production is about ±10%. 

In order to evaluate what neutrino energies are relevant 
for the T2K v e appearance search, one has to look at the 
typical energies of pions that constitute background. The 
oscillation signal from — » v e is peaked around 0.7 GeV 



of the the neutrino energy 14| . This means that danger- 
ous 7r°s have energy in this region. Considering energies 
between 0.4 and 0.9 GeV and using Nuance simulation 
to translate them into incident neutrino energies, we get 
neutrino energy range 1 to 4 GeV. We can then conclude 
that neutrino of such energies are the most important 
in the context of v e appearance searches in the T2K ex- 
periment. Therefore the expected uncertainty varies be- 
tween 6% and 10%. However, it is important to note 
that the errors discussed here come only from uncertain- 
ties in form-factors and do not take into account other 
model approximations, e.g. nonresonant background in 
pion production or an influence of matter on the width 
of A resonance. 

We also notice that Nuance predicts higher (with re- 
spect to NuWro) SPP cross-sections on nuclear targets. 
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FIG. 7: cross-sections of SPP production on water. From 
top to bottom - CC ln + production absolute cross-section, 
NC 7T° production absolute cross-section, CC ir + production 
cross-section normalized to Nuance results without NEO, and 
NC 7T° production cross-section normalized to Nuance results 
without NEO. NuWro points show la error contours. Nuance 
points are shown here for reference. 
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FIG. 8: The ratio CC17T+/CCQE on the CH 2 target. Data 
points are taken from fliil ]. The dipole fit of this paper was 
used in the NuWro Monte Carlo event generator with QE 
axial mass M% E = 1.03 GeV. 



This is probably caused by different assumptions for pion 
absorption cross-section at higher values of kinetic en- 
ergy. This is a region where the experimental data for 
pion nucleus absorption are missing and Monte Carlo pre- 
dictions have to rely on theoretical assumptions. 



7.2. Comparison with measured CC17T+/CCQE 
ratio 



Recently there has been much controversy about the 
overall normalization of the neutrino cross-section. The 
MiniBooNE and NOMAD experiments report very differ- 
ent values of the QE axial mass (Ma = 1-35 ± 0.17 GeV 
k = 1.07±0.07 or M A = 1.37±0.12 GeV with k = 1.0 and 
Ma = 1-05 ± 0.08 GeV respectively). For this reason the 
most reliable measurements are free from normalization 
uncertainties, and an interesting example of such mea- 
surements is the recent CC17T+/CCQE ratio from Mini- 
BooNE. In the Ref. [l6| two sets of data are presented, 
one for lir + on the nuclear target CH2 and the second 
with evaluation on the same ratio on isoscalar nuclear tar- 
get. We compare with the first set of data which is more 
independent on models contains in MiniBooNE's Monte 
Carlo generator of events. The value of the QE axial mass 
used in the simulation was Ml E = 1.03 GeV. In Fig. Q 
we see that we obtained quite good agreement with the 
MiniBooNE data. The MiniBooNE CC 1tt+ /CCQE data 
imply that the values of QE and A production the ax- 
ial masses are correlated. If one increases in NuWro the 
value of the axial mass by 30% as suggested by the Mini- 
BooNE QE data, then in order to keep agreement with 
the measured ratio, the value of axial mass in C^(Q 2 ) 
would have to be accordingly increased. 
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8. FINAL REMARKS 

We have demonstrated that old bubble chamber ANL 
and BNL SPP data are self-consistent and restrictive 
enough to enable one to extract quite precise information 
about the axial form-factor C^(Q 2 ). An important ingre- 
dient of our reanalysis was the inclusion of normalization 
(flux) uncertainty. We also took into account deuteron 
nuclear effects but their impact on the final result was 
surprisingly small due to the interplay with normaliza- 
tion factors. 



It can be shown that \ is distributed with: 

D 



The evaluation of nuclear effects presented in Ref. [12 1 
was based on nonrelativistic approximation and it would 
be interesting to perform exact computations and inves- 
tigate how much this would modify the numerical results 
of this paper. 

It would also be interesting to try to extract infor- 
mation about the nonresonant background contained 
in ANL and BNL SPP data on bound neutron with 
deuteron effects included. Here the situation becomes 
more complicated because deuteron FSI effects like Pauli 
blocking can play an important role. 
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APPENDIX A: PARAMETER GOODNESS OF 
FIT 

Suppose that D independent data sets are analyzed. 
Every data set contains N r (r = 1,2,..., D) bins (for a 
review, see Ref. [ill). Let P be the total number of 
parameters which are fitted. One can construct the \r 
for the i-th data set Thus the total \ 2 for the global fit 
reads 



D 

2 — \ ~" v v 2 

Xiot / j Xr 

r=l 



(Al) 



The idea of the parameter goodness of fit is to consider 
a redefined x 2 , namely 



D 

— 2 2 2 

X = Xtot - l^Xr .nun " 

r=l 



NDF C = Pr — P- 



(A3) 



degrees of freedom. 

Then, the parameter goodness of fit is defined as: 

PGoF = CL(x,NDF c ). 



(A4) 



Notice that \ has a minimum at the same point as Xto 
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